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We study the normal state electronic excitations probed by angle resolved photoemission spectroscopy 
(ARPES) in Bii.6Pbo.4,Sr2Cu06 (Bi2201) and Bi2Sr2CaCu20s+s (Bi2212). Our main goal is to establish 
explicit criteria for determining the Fermi surface from ARPES data on strongly interacting systems where 
sharply defined quasiparticles do not exist and the dispersion is very weak in parts of the Brillouin zone. 
Additional complications arise from strong matrix element variations within the zone. We present detailed 
results as a function of incident photon energy, and show simple experimental tests to distinguish between 
an intensity drop due to matrix element effects and spectral weight loss due to a Fermi crossing. We reiterate 
the use of polarization selection rules in disentangling the effect of umklapps due to the BiO superlattice in 
Bi22i2. We conclude that, despite all the complications, the Fermi surface can be determined unambiguously: 
it is a single large hole barrel centered about (vr, tv) in both materials. 

PACS numbers: 71.25.Hc, 74.25.Jb, 74.72.Hs, 79.60.Bm 
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I. INTRODUCTION 

The electronic structure and Fermi surfaces of con- 
ventional metals have been studied in great detail by 
ARPES iL The question of the determination of the Fermi 
surfacem3 by ARPES in the normal state of the high 
Tc superconductors is of great interest, especially since 
other Fermi surface probes (like de Haas- van Alphen and 
positrons) have not yet yielded useful information on the 
Fermi surface of the planar Cu-0 states. However, this 
question is not a trivial one, since these materials are 
strongly correlated, and likely not Fermi liquidSp,exhibit- 
ing very broad, ill-defined electronic excitationaJo. 

The determination of the Fermi surface by ARPES in 
these systems is further complicated by the very small 
dispersionJn-the vicinity of the (tt, 0) point of the Bril- 
louin zonm'Lm, and by strongly k-dependent photoemis- 
sion matrix elements which lead to intensity variations 
which have nothing to do with Fermi crossings. Bi2212 
has an additional complication: final state diffraction 
of photoelcctrons by the Q =riWr^'''i 0.2l7r) structural 
modulation in the BiO layera3l3'E2l. The combination 
of all these effects, if not treated correctly, can be a 
source of confusinn and lead to apparently contradictory 
conclusionsoild lij, even though the data between various 
groups are completely consistent with one another. 

It is therefore important to establish the criteria for un- 



ambiguously extracting the Fermi surface from ARPES 
data. This is the main goal of the work presented here, 
in which we study near-optimal and overdoped samples 
oi Bi2Sr2CaCu20s+s (Bi2212) and Bii,ePbo,4Sr2CuOe 
(Bi2201) using a range of incident photon energies from 
17-60 eV. Our extensive studiz-Jeads us to the same con- 
clusion as our previous workO'Ej, namely that the Fermi 
surface consists of a single hole barrel centered around 
(tt, tt), the most antibonding point in the Brillouin zone. 
Furthermore, the Fermi surface is consistent with the 
Luttinger count with its volume scaling as one plus the 
number of doped holes. 

The paper is organized as follows. In Section II, we 
describe the samples and experimental details. We begin 
with ARPES data over a wide (6eV) scale describing the 
full valence band in Section HI, and then turn to low en- 
ergy, near Ep features in the rest of the paper. Section IV 
contains a brief description of polarization selection rules 
and their experimental implications. We then discuss 
in some detail criteria for determining the Fermi surface 
in Section V. In Section VI we present data on Bi2212 
and Bi2201 and show how the various criteria proposed 
in Section V fare in determining Fermi crossings. We 
find that the symmetrization method for inferring when 
the spectral function peak goes through the chemical po- 
tential is a very powerful tool and works even when the 
energy distribution curves (EDCs) are broad and weakly 
dispersive. We discuss in Section VII and VIII the useful- 



ness and limitations of using the integrated intensity to 
determine k-p , emphasizing the importance of the photon 
energy and k-dependence of the ARPES matrix elements. 
By analyzing data obtained at different incident photon 
energies, we explicitly show how one can experimentally 
separate matrix element effects from those due to changes 
in the momentum distribution. In Section IX, we finally 
turn to the Fermi surface in Bi2212 where, in addition to 
all the issues discussed above for Bi2201, one also needs 
to be careful about BiO superlattice effects. Polarization 
selection rules are exploited to disentangle superlattice 
effects from the intrinsic Cu02 electronic structure. We 
conclude in Section X. 

An Appendix contains some further technical details 
related to the symmetrization procedure. 



II. EXPERIMENTAL DETAILS 

Our experiments are on very high quality single crys- 
tals of Bi2212 and Bi2201 grown by the traveling solvent 
floating zone method with an infrared mirror furnace, 
with low defect densities, as it can be appreciated from 
the high resolution x-ray diffraction rocking curve shown 
in Fig. 1. The samples are labeled by their doping levels 
(OPT for optimal doped and OD for overdoped) together 
with their onset T^.. 

The as-grown Bi2212 samples are slightly overdoped 
and have Tc = 87K with a transition width of IK as 
determined by a SQUID magnetometer. These samples 
are most stable in terms of their photoemission charac- 
teristics. We look at Bi2201 samples in the doping range 
from OD23K to heavily ODOK. The samples are cleaved 
in situ, and have optically flat surfaces as measured by 
specular laser reflections. It is absolutely essential to 
characterize the flatness of the surface on which ARPES 
experiments are done. Another measure of the sample 
qualitiL within ARPES, is the observation of "umklapp" 
bandstl in the electronic structure of Bi2212 samples, due 
to the presence of a structural superlattice modulation in 
the BiO layer. Since the structural superlattice has a pe- 
riodicity of ^ 5 times the unit cell, very good long range 
order is required for its observation. 

The experiments were performed at the Synchrotron 
Radiation Center, Wisconsin, using a high resolution 4-m 
normal incidence monochromator with a resolving power 
of 10"' at 10^^ photons/s. The samples are carefully ori- 
ented in the sample holder to an accuracy of 1° by Laue 
diffraction, and the orientation is further conflrmed by 
the observed symmetry of sharp ARPES features around 
high symmetry points, as described below. 

Some of the data analyzed below (in particular Figures 
5, 11, 14 and 16) were obtained using a Scienta analyzer, 
at a variety of incident photon energies, with an energy 
resolution of 16 meV and a high k-resolution better than 
0.01 A~^ at 22 eV photon energy. The detailed analysis 
presented in this paper leads to conclusions which are 



completely consistent with the recent high k-resolution 
results of our groupEj (emphasizing low tempeifature data 
on Bi2212), as well as that of Borisenko et al.ll3. 

For the Brillouin zone of Bi2212 and Bi2201, we use 
a square lattice notation with TM along the CuO bond 
direction, as shown in the insets of Fig. 2. F = (0,0), 
M ~ (tt, 0), X = (tt, — tt) and Y = (tt, tt) in units of 
1/a*, where a* = 3. 83 A is the separation between near 
neighbor Cu ions. (The orthorhombic a axis is along X 
and b axis along Y). 

An example of how ARPES is used in sample align- 
ment is shown in Fig. 2, where spectra are shown along 
the (0, 0) — (27r, 0) and the (tt, tt) — (tt, — tt) high symme- 
try lines of the Brillouin zone. This symmetry is reflected 
in the position of the peak in the spectra in Fig. 2, and 
allows us to accurately flnd the surface normal, and the 
angle of the sample about this normal, completely deter- 
mining the momentum k of the outgoing electron. Note 
that this alignment procedure only makes, use of the sym- 
metry properties of the peak positionalJ, and does not 
require a knowledge of the Fermi surface (indicated by 
the curves in the top panels of Fig. 2). 



III. THE VALENCE BAND 

Our main focus will be on near Ep electronic structure, 
but we begin with a brief discussion of angle-resolved 
photoemission from the complete valence band of Bi2212. 
This covers an energy range of approximately 6 eV from 
the Fermi energy (the small peak near zero) to the bot- 
tom of the valence band. The electronic structure can 
be divided into three groups, as indicated in Fig. 3a: the 
most bonding Cu02 state is at the bottom of the valence 
band (the peak at 6 eV), the antibonding state is at the 
Fermi energy, and the non-bonding states are in between. 
The "lump" in the middle also includes states from ele- 
ments in the structure other than Cu and O. But since 
these layers are insulating, the corresponding states do 
not cross the Fermi energy. 

By varying the in-plane momentum k, one can map 
the complete electronic structure of the valence band, as 
shown in Fig. 3a. These curves were obtained without flt- 
ting peaks to the data. Instead, the second derivative of 
the observed spectra, as shown in Fig. 3b, was taJsen and 
plotted as a grey scale without any modificationstj. Two 
considerations apply: the energy step in these spectra is 
only 30 meV, and therefore the details of the dispersion 
near the Fermi energy are not clear, and the spectra were 
obtained with a particular polarization of the photons, so 
that not all states show optimal intensity. Nonetheless, 
the most noteworthy features are the most bonding and 
antibonding states, highlighted by thick dark lines. In 
the remainder of the paper, we will focus exclusively on 
the antibonding states in the region near the Fermi en- 
ergy. 



IV. SELECTION RULES 



V. FERMI SURFACE CRITERIA 



We now focus on the states crossing the Fermi energy 
in Fig. 3a and show how one can determine the symmetry 
of the initial state in ARPES. The ARPES intensity is 
governed by the (square of the) dipole matrix element 
Mfi connecting the initial state |^i) to the final state 
IV'/), given by |(^/| A • p |^i)| , where A is the vector 
potential of the (linearly polarized) incident photon, and 
p is the momentum operator. 

We use a siroDle version of the selection rules proposed 
by HermansonEJ. Let the photon beam be incident along 
a plane of mirror symmetry of the sample {M). If the 
detector is placed in the same mirror plane as shown in 
Fig. 4a, then the final state ^pf must be even with re- 
spect to reflection in A1 , because if it were odd the wave 
function would vanish at the detector. The dipole transi- 
tion is allowed if the entire matrix element has,an overall 
even symmetry. Thus two possibilities arisecJ. First, if 
the initial state ipi is even with respect to A4, then the 
light polarization A must also be even, i.e. parallel to 
Ai. Second, if the initial state is odd with respect to Ai, 
then A must also be odd, i.e. perpendicular to A4. This 
can be summarized as: 



(V'/|A-p|V'») 



Ipi even 
-tpi odd 



+ 1+) 



A even 
A odd 



(1) 



Consider hybridized Cuid — 02p initial states, as 
shown in Fig. 4b, which have a dx^_y2 symmetry about 
a Cu site. These states are even with respect to (0, 0) — 
(tt, 0) (i.e. the plane defined by this symmetry axis and 
the z-axis) and odd with respect to (0, 0) — (tt, tt). There- 
fore, measurement along the (0, 0) — (tt, 0) direction will 
be dipole-allowed (forbidden) if the polarization vector 
A is parallel (perpendicular) to this axis. Fig. 4c shows 
that, consistent with an initial state which is even about 
(0,0) — (tt, 0), the signal is maximized when A lies in 
the mirror plane and minimized when A is perpendicu- 
lar to this plane. (The reasons for non-zero intensity in 
the dipole forbidden geometry are the small, but finite, 
k- window of the experiment and the possibility of a small 
misalignment of the sample) . Similarly, we have checked 
experimentally that (for Bi2212 in the Y-quadrant where 
there are no superlattice complications) the initial state 
is consistent with odd symmetry about (0, 0) — (tt, tt). 

While the dipole matrix elements are strongly pho- 
ton energy dependent, the selection rules are, of course, 
independent of photon energy. This has been checked 
by measurements at 22 eV and 34 eV. All of these re- 
sults are consistent with the fact that we are probing 
Cu3d — 02p initial states with dx2-y2 symmetry. In ad- 
dition, as we shall emphasize below, the selection rules 
can be exploited to one's great advantage in disentan- 
gling the main Cu02 "band" from its umklapp images 
due to the superlattice in Bi2212. 



Many criteria have been used for determining the 
Fermi surface in the past without a clear discussion of 
the conditions under which they are applicable. We will 
present three criteria here: (A) one based on dispersion of 
the EDC spectral peaks through the chemical potential, 
(B) a second one based on the peak of the spectral func- 
tion inferred from symmetrized data, and (C) a third one 
based on rapid changes in the momentum distribution. 
In the following Sections we will show how these criteria 
fare when applied to experimental data. Other criteria, 
not discussed in this paper, will be briefly alluded to at 
the end of this Section. j_. 

The ARPES intensity is given byB 



I{k,uj) = Io{k;,y;A)f{uj)A{-k,uj) 



(2) 



for a quasi-two-dimensional system, assuming validity of 
the impulse approximation. Here k is the in-plane mo- 
mentum, Lu is the energy of the initial state measured rel- 
ative to the chemical potential, /(w) = l/[exp((jj/r) -I- 1] 
is the Fermi function, and the one-particle spectral func- 
tion A(k,w) = (-l/7r)9TOG(k,w-FJ0+). The prefactor /o 
is proportional to the dipole matrix element \Mfi\ and 
thus a function of k and of the incident photon energy 
hv and polarization A. It is also important to remember 
that the experimentally observed EDC involves a convo- 
lution of the intensity of Eq. (0) with the energy resolu- 
tion function and a sum over the momentum resolution 
window. There is also an additive (extrinsic) background 
contribution to the EDC, however this has little effect on 
Fermi surface determination since, the background is neg- 
ligible at the chemical potentiate. 

We first discuss the simplest case of non-interacting 
electrons which have infinitely sharp energy levels leading 
to a spectral function A(k, w) = S{uj — e^)- A Fermi 
surface crossing kj? is then defined by the location in 
k-space where the sharp peak of the spectral function 
crosses the chemical potential {u; — 0), i.e., ekj? — 0- It 
is also useful to look at the momentum distribution 



ri(k) 



dujf{uj)A{'k,uj). 



(3) 



For non-interacting electrons n(k) = /(^k): the Fermi 
function. At T = then, the momentum distribution 
shows a jump discontinuity at kp. At low temperatures 
there is no singularity in n(k) but only a rapid variation 
in the vicinity of kp. 

The case of interacting electrons at finite temperatures 
is much more interesting. The energy levels are now 
broadened and shifted by the self-energy S = S' -I- iS", 
with G~^(k, w) = w — Ek — S(k, w). Thus the spectral 
function is given by 



1 



Aik,u;) = - 



|E"(k,u;) 



^[c^-ek-E'(k,o.)]^ + [S"(k, 



(4) 



The electronic dispersion is now given by tracking the 
peak of the spectral function. We define a Fermi surface 
crossing by the k-point at which the spectral function 
peak crosses the chemical potential {uj = 0) 



sReG-i(kF,w = 0) = 0. 



(5) 



This definition agrees with the standard definition of 
the Fermi surface at T = in an interacting Fermi sys- 
tem wiiich can be described by Landau's Fermi liquid 
theorynJ. For this case, there is an additional (equiva- 
lent) characterization of the Fermi surface in terms of a 
discontinuity in 7i(k) at T = 0. 

However, it should be stressed that the discussion 
above in terms of the spectral function is very general, 
and not necessarily limited to Fermi liquids. We propose 
to use the above definition in terms of the peak of the 
spectral function, which is valid even at finite tempera- 
tures, and use it to define the Fermi surface for high Tc 
superconductors even though the spectral peaks above 
Tc are too bcoad for the system to qualify as a Landau 
Fermi liquidQ. The question of whether the T — mo- 
mentum distribution shows any singularity or not cannot 
be addressed experimentally, since at T = one is not in 
a normal state but rather a broken symmetry state. We 
do find, however, that the Fermi surface we experimen- 
tally determine above Tc encloses a number xif electrons 
which are consistent with the Luttinger counted of (H-x), 
where x is the hole doping. 

Let us now discuss in detail how these characterizations 
of the Fermi surface will be used in practice to determine 
ki;-. The first method (A) is simply to look at the disper- 
sion of the peaks of the measured EDCs and determine 
from this when the peak position crosses the chemical 
potential. There are two caveats to this method. First, 
the peak of the EDC does not in general coincide with 
the peak of the spectral function A(k,uj). As can be 
seen from Eq. (||), if there is a broad spectral function 
A centered about a; = 0, then the peak of the EDC will 
be at w < 0, produced by the Fermi function chopping 
off the peak of A, in addition to resolution effects. This 
can readily be seen in the data, as will be discussed in 
connection with Fig. 5(c) below, and can be corrected 
for under favorable circumstances. The second problem 
with method (A) is that it may be difficult to use in 
cases where the dispersion is very weak, as for instance 
near the (tt, 0) point in the cuprates. We should note 
that with sufficiently fine k-sampling, these problems are 
minimized, as shown recently by us in ref. n^. However, 
as we show here, even in the absence of such data it is 
possible to make progress. 

We turn to the symmetrization method (B) which al- 
lows us to overcome both the limitations of the meiJaod 
(A). This method was originally introduced by ua£3 as 
a means of "dividing out the Fermi function" from the 
EDC and directly infer the spectral function A. In order 
to determine A(ki?,a;) one had to assume particle-hole 
symmetry on a low energy scale. However, we now use 



this idea for a different purpose, namely Fermi surface 
determination. As we show below we do not need any 
assumptions about p-h symmetry to determine hp- 

For an arbitrary k, we define the symmetrized ARPES 
intensity by 

Isym{k,Uj)^I{k,Uj)+Iik,-ij) (6) 

For simplicity, wc will ignore resolution effects here, and 
we use eq. (||) on the right hand side of eq. (^). (The 
effect of energy and momentum resolution convolutions 
are discussed in the Appendix.) Our goal is to use /gym 
to determine kp at which A has a maximum at uj = 0. 

Analyzing the symmetrized intensity about lu — 0, for 
any k, we see that dlsym/ dLo{uj 



d^h. 



duj^ 



j=0 



'sym 

fA 



j=0 



0) = and 

I dA 

T dio 



(7) 

For k 7^ kp, the second term on the right hand side 
dominates at sufficiently low temperature. For an oc- 
cupied state, it is easy to see that dA/duj{LL! = 0) < 0, 
so that d^/sym/dw^(w = 0) is positive. Thus the sym- 
metrized intensity will exhibit a local minimum, or a dip, 
at ll!=0 for an occupied k-state. Conversely, for an unoc- 
cupied state dA/duj{uj = 0) > 0, and the left hand side 
of eq. (0) is negative, leading to a local maximum in the 
symmetrized intensity at lu=0. Precisely at k = kp, the 
spectral function has a maximum at w = 0. Thus the 
second term of eq. (0) vanishes and the first term leads 



to d Isym/duj {lu 
mum in Igym at '^ 



0) < yielding a peak or local maxi- 

0. 
In practice, symmetrization is used to determine the 
Fermi crossing kp as follows. We symmetrize all EDCs 
along a cut in k-space and identify kp as the boundary in 
momentum space between where symmetrized data have 
a dip (local minimum) versus a peak (local maximum) at 
UJ — 0. This will be demonstrated in detail below, where 
the symmetrization estimate for kp from ARPES data is 
also compared with other estimates, wherever possible, 
and found to agree. 

Before turning to the data (in the next Section), we 
show one example of a simulation which illustrates sym- 
metrization with resolution effects included. In Fig. 16 
we plot symmetrized intensities for five k points along 
a certain cut in k-space. We see that for occupied 
states, corresponding to the first two k points, the spec- 
tral functions peak at a; < 0, and thus the correspond- 
ing symmetrized spectra show dips, or local minima, at 
UJ = 0. As kp is approached, this minimum gets shal- 
lower. At kp, half the peak is being chopped off in the 
EDC. Therefore, upon symmetrization it is completely 
restored. (Note we did not build in any matrix element 
effects in this simulation). Once you go beyond kp, as in 
the last two curves, more than half the peak is chopped 
off in the EDC, so there is an intensity drop in the sym- 
metrized plot. In addition, one observes that unoccupied 
states continue to exhibit a peak at w = in their sym- 
metrized intensity in Fig. 16, as can be understood from 
the preceding discussion. 



We note that symmetrization makes precise, or formal- 
izes, a rough criterion often used by ARPES practition- 
ers: at kp a vertical line through the chemical potential 
intersects the midpoint (half-the maximum intensity) of 
the leading edge of the EDCu. It should be emphasized 
a very precise determination of the chemical potential 
(to = 0) is necessary to determine hp via symmetriza- 
tion. 

pEipally we turn to another method (C) devised by 
ufO'ESl, based on the sum rule Eq. (g) relating the energy- 
integrated ARPES intensity to the momentum distribu- 
tion n(k). In principle, the rapid variation of n(k) offers 
a very direct probe of the Fermi surface, which is again 
not restricted to Fermi liquids. (The T = momentum 
distribution for known non- Fermi liquid systems, such as 
Luttinger liquids in one dimension, shows an inflection at 
ki?.) In several cases we have demonstcated the useful- 
ness of this method in our earlier workua where \i.p was 
estimated from the location of max |Vk'^(k)|. The same 
methocLiias also been successfully used later by other 
authorsBB. 

However, there is an important caveat to keep in mind: 
one does not measure f^(k), but rather the integrated 
intensity 



/(k) 



+ 00 



dujl{k,uj) = /o(k;i/;A)n(k). 



(8) 



Potential problems for Fermi surface determination 
can, and do, arise from the k-dependence of prefactor 
/o(k;t/;A) due to ARPES matrix elements. In Sections 
VII and VIII below we discuss in detail how to distin- 
guish k-dependences of the integrated intensity coming 
from n(k) and from the matrix elements. 

To conclude our discussion of Fermi surface criteria, 
we note that we will restrict ourselves here to the normal 
state. We will not dis cu s s jn this paper the notion of the 
"minimum gap locus"c30 in a gapped state (either be- 
low Tc or in the pseudogap regime) which is a measure- 
ment of the underlying Fermi surface that got gapped 
out. We also mention, for completeness, two other meth- 
ods of Fermi surface determination, which we will not 
discuss in thisjiaper. The first one exploits the approxi- 
mate sum ruleEJ that dn{\ip) / dT = 0, i.e., the integrated 
intensity at ki? is independent of temperature. The sec- 
ond method uses the constant uj scan, or momentum dis- 
tribution curve (MDC), as a function of k, at a; = 0; see 
ref. 51 



VI. FERMI CROSSING FROM 
SYMMETRIZATION 

We will show below that the symmetrization method 
provides a simple and general way of determining a Fermi 
crossing, even when the dispersion is very small. As dis- 
cussed above, we will identify k^? as that k for which the 



symmetrized data first shows a clear peak at the chemical 
potential (w = 0). 

We begin with the simplest Fermi crossing along the 
zone diagonal, (0, 0) — > (tt, tt), where the electronic struc- 
ture shows rapid dispersion. The data for a Bi2212- 
OD88K sample shown in Fig. 5(a) were obtained with 
an incident photon energy of 22 eV and polarization A 
parallel to the (0, 0) - (tt, -tt) axis. These T = 180K 
spectra are extremely broad^nd not at all consistent 
with a Fermi liquid pictureaE3. We first determine the 
Fermi crossing \<ip using the dispersion of the observed 
peak. In Fig. 5(c) we plot the observed peak positions of 
the data in (a) . One can see that in the narrow k interval 
plotted in Fig. 5(c), the dispersion is linear over much of 
the range, deviating from linearity due to the effect of 
the Fermi function (as discussed above). The extrapo- 
lation of the linear part crosses the chemical potential 
at k = (0.3757r, 0.3757r), which is then the estimated kj? 
from the dispersion. 

In Fig. 5(b) we plot the symmetrized intensities 
/sym(k,ti') obtained from the data of Fig. 5(a). From 
this we sec that the inferred spectral function is peaked 
at w = also at k = (0.377r,0.377r). We thus find that 
in this case the dispersion and symmetrization methods 
give identical \ip estimates. 

In the following we will first look at the simpler case of 
Bi2201. In the Pb-doped Bi2201 compounds there are no 
observable complications arising from umklapp bands as 
in oxygen-doped Bi2201 or Bi2212 (which is discussed in 
more detail in Section IX) . We now move along the Fermi 
surface in Fig. 6, where the plotted data were obtained 
for an OD23K sample with an incident photon energy of 
22 eV and polarization A parallel to the (0, 0) — (tt, 0) 
axis. One can observe a clear trend: as the (tt, 0) — (tt, tt) 
line is approached, the dispersion becomes very small. In 
the vicinity of the (tt, O) it becomes very difficult to use 
the dispersion criterion to determine kj?. Nevertheless, 
symmetrized data provide completely unambiguous re- 
sults as can be seen from Fig. 7. From the top panels 
(a,b, and c) of Fig. 7 one can determine precisely the 
Fermi crossing from the k-point at which the spectral 
function inferred from symmetrized data first peaks at 
the chemical potential. 

Two points should be noted about the k^? estimate 
from symmetrization. First, just before approaching \<ip 
from the occupied side, we expect that the symmetrized 
data will show two peaks and a small dip, which would 
be broadened by resolution into a rather flat topped sym- 
metrized spectrum. Second, there should be an intensity 
drop upon crossing \<ip in the symmetrized spectrum, as- 
suming that the matrix elements are not strong functions 
of k. Both of these effects are clearly seen in the data. 

It is equally important to be able to determine the 
absence of a Fermi crossing in a cut in k-space as shown 
in Fig. 7(d,e). In this respect, the raw data along (0, 0) — 
(tt, 0) — (27r, 0) is quite difficult to interpret, since the 
"band" flattens while approaching (7r,0), and remains 
extremely close to Ep. Nevertheless, it is simple to see 



that at no point along this cut do the symmetrized data 
show a peak centered at Ef, thus estabhshing the absence 
of a Fermi crossing along this cut. 

The important conclusion from this discussion is that 
for Pb-doped Bi2201 one can continuously follow a Fermi 
surface contour which traces a hole barrel centered at 
(tt, tt). Analysis of a large set of data using the sym- 
metrization method allows an unambiguous determina- 
tion of the Fermi surface crossing with high accuracy, 
even in the unfavorable case of broad peaks with weak 
dispersion. 



VII. MATRIX ELEMENTS 

As discussed at the end of Section V, great care must 
be used in determining a Fermi crossing from the in- 
tegrated intensity which is the momentum distribution 
n(k) multiplied by the prcfactor /o(k;i^;A). A loss of 
integrated intensity as a function of k can arise either 
from a drop in n(k) related to a Fermi crossing, or from 
the k-dependence of the matrix elements in Iq . 

One possibility is to have a priori information 
about the |fflatrix elements from electronic structure 
calculation^. But as we now show, even in the absence 
of such information, one can experimentally separate the 
effects of a strong k- variation of the matrix element from 
a true Fermi surface crossing. The basic idea is to ex- 
ploit the fact that changing the incident photon energy 
one only changes the ARPES matrix elements and not 
the momentum distribution of the initial states. 

In Fig. 8(a) we shows cuts along (0,0) — (tt, 0) and 
(tt, 0) — (tt, tt) obtained at a photon energy of 34 eV (to 
be contrasted with the data in the previous two figures 
for the same sample at 22 eV). The symmetrized data 
are shown in Figs. 8(b) and (c) from which we see results 
entirely consistent with those obtained at 22 eV. There 
is no Fermi crossing along (0, 0) — (tt, 0) since the sym- 
metrized data in Fig. 8(b) never show a peak at a; = 0. 
Turning to the (tt, 0) — (tt, tt) direction, the symmetrized 
data (Fig. 8(c)) do show a Fermi crossing occurring at 
k=0.127r, in agreement with the data obtained at 22 eV 
(k=0.147r). 

However, there is an important difference between the 
data sets at 22 eV and 34 eV photon energies, which 
can be appreciated in Fig. 9, where the integrated inten- 
sity along the two directions is displayed. While at 22 
eV the maximum intensity occurs close to (tt, 0) and de- 
creases both towards (0,0) and (7r,7r), the data taken at 
34 eV show a strong depression of intensity on approach- 
ing (7r,0), resulting in a shift of the intensity maximum 
away from (tt, 0). This loss of intensity cannot be in- 
terpreted as a Fermi crossing, since at no point in the 
symmetrized data from (0,0) — > (7r,0) (Fig. 8) is there 
a peak centered at w = 0. (In fact, this loss of intensity 
close to (tt, 0) is responsible for the reduced signal-to- 
noise ratio in the 34 eV data in Fig. 8 compared with the 



22 eV data in Figs. 6 and 7). 

We would like to attribute this loss in intensity around 
(tt, 0) at 34 eV, and in fact the entire variation seen in 
Fig. 9(a), to strong k-depcndcnt matrix element effects. 
A direct proof is found from the data: the EDCs at the 
same point in the Brillouin zone obtained at the two 
different photon energies exhibit exactly the same line- 
shape, i.e. one can be rescaled onto the other as shown 
in Fig. 9b. 

We emphasize that the results of Fig. 9 imply that the 
photon energy dependence of the ARPES data is not a k^ 
dispersion effect. If this were the case, different incident 
photon energies would be probing initial states with dif- 
ferent kz values. However, the scaling of Fig. 9(b) proves 
that it is the same two dimensional (fc^-independent) ini- 
tial state which is being probed, and the photon energy 
dependence arises from the different final states that the 
matrix element couples to. 

To further illustrate the role of matrix elements, we 
present a model calculation and compare it with the data 
of Fig. 9. The purpose of this exercise is to determine 
whether or not the data are consistent with a matrix 
element variation with k. In Fig. 10(a) we use dotted 
lines to show the dispersion of a model ARPES spec- 
trum along (0,0) to (7r,0) where the "band" approaches 
Ep near (tt, 0) without a Fermi crossing. (The dispersion 
is chosen from the tight binding fit of Ref. p3 to ARPES 
data on Bi2212 for illustrative purposes even though we 
will compare it to data on Bi2201). The full curves in 
Fig. 10(a) show the effect of the matrix element variation 
on the model spectra, simulated by /o(k) = sin (O.Gfc^,) 
along (0,0) to (7r,0). This is a simple phenomenologi- 
cal matrix element (squared) which satisfies the follow- 
ing properties: it vanishes at the F point, as dictated by 
symmetry, and then has non-monotonic behavior along 
F-M with a peak away from the M point, similar jto that 
obtained by detailed band-structure calculationsEfl. We 
emphasize that, beyond this, no deep meaning should be 
attached to the simple analytical form used. 

Fig. 10(b) shows the momentum dependence of the in- 
tensity integrated over the (large) energy range of -350 
to +50 meV as a solid line, and over the (narrow) en- 
ergy range -50 to -1-50 meV as a dashed line. One can 
see that this simple example shows remarkable agreement 
with the measured intensity using the same integration 
ranges, shown in Fig. 10(c). These results are easy to un- 
derstand. The data with the large integration range first 
increase in intensity simply following the matrix element 
variation. The data with integration over a small energy 
range however shows a different intensity behavior. They 
start to increase for k-values higher than 0.657r because 
only when the peak is closer to the Fermi energy does 
it contribute to the integrated intensity. It then again 
decreases rapidly because of the strong decrease of the 
matrix element, as in the previous case. 

However, not recognizing the role of matrix elements, 
some authors have ascribed the differences between the 
22 eV and 34 eV photon energy data to additional "mys- 



terious" states around (tt, 0)li3. We feel that there is no 
necessity to invoke such states or to assert changes in 
Fermi surface topology with photon energy. In conclu- 
sion, all our data on Bi2201 when analyzed using the 
methods described above indicate a large Fermi surface 
centered around (tt, tt) independent of the photon energy. 



VIII. THE MOMENTUM DISTRIBUTION 

Despite the matrix element issues discussed above, it is 
nevertheless interesting to study the integrated intensity 
(Eq. (H)) for a dense data set in the entire zone. The re- 
sults obtained are shown in Fig. 11 for an ODOK sample. 
In the top panels Figs. 11(a) and (c) we show the inte- 
grated intensity J(k) around the (tt, 0) point obtained at 
two different photon energies: 22 eV and 28 eV respec- 
tively. In the lower panels Figs 11(b) and (d), we plot the 
magnitude of the logarithmic gradient: |Vk/(k)|//(k) 
which emphasizes the rapid changes in the integrated 
intensity. The logarithmic gradient filters out the less 
abrupt changes in the matrix elements and helps to fo- 
cus on the intrinsic variations in n(k). This can be seen 
from the fact that the integrated intensities in the top 
panels are quite different for 22 eV and 28 eV, while the 
logarithmic gradients are much more similar. 

The Fermi surface can be clearly seen as two high in- 
tensity arcs curving away from the (7r,0) point. Modulo 
matrix element effects, the results obtained at the two 
different photon energies are quite similar. Moreover, 
the Fermi surface estimated by this method is in good 
agreement with the one obtained from the symmetriza- 
tion analysis above. 



IX. THE FERMI SURFACE OF Bi2212 

We now turn to Bi2212 where, in addition to all the 
issues discussed above for Bi2201, there is an added com- 
plication due to the presence of umklapp bands aris- 
ing from the superlattice modulation with wavevector 
Q = (0.2l7r,0.2l7r) in the BiO layers. Another differ- 
ence with Bi2201, which we will not address here, is that 
Bi2212 has a Cu02 bilayer, we only mention that no bi- 
layer splitting is observed in the ARPES data on Bi2212 
as discussed in detail in Ref. 1^. In this Section we will 
first review the effect of the superlattice on the electronic 
structure probed by ARPES, emphasizing the usefuJness 
of polarization selection rules. In an earlier lettem, we 
had shown data along the principal axes of Bi2212. Here 
we present cuts throughout the Brillouin zone, analyzed 
using the symmetrization and integrated intensity meth- 
ods discussed above, together with a detailed study of 
the photon energy dependence of the matrix elements. 
All of these new data and their analysis substantially 
strengthen our earlier conclusions. (See also ref. |l5[). 



In Fig. 12 we show the electronic structure from EDC 
peaks (in panel (c)) and Fermi surface crossings (in panel 
(b)) determined from data at incident ^photon energies of 
19 and 22 eV for an OD87K samplel. The dark lines 
in the bottom panel of Fig. 12 are a fit of the intrinsic 
planar Cu02 electronic structure, ek, which we call the 
"main band"; see Ref. p|j3^ for details. The lighter lines 
are simply obtained by plotting ek±Q, where Q is the 
superlattice vector, and it is very important to note that 
these lines provide an excellent description of the data 
points tlw,t do not lie on main band. The data strongly 
suggestES that these additional "umklapp bands" arise 
due to diffraction of the outgoing photoelectron through 
the BiO superlattice, which leads to "ghost" images of 
the electronic structure at ek±Q- 

From the point of view of the present discussion of the 
Fermi crossings, it is very important to establish con- 
clusively that the crossings U4 and U5 along (0, 0) to 
(tt, 0) shown in the middle panel of Fig. 12 correspond 
to umklapp "ghost" images and not to the "main band" . 
The case of the Fermi crossing U5, closer to (0,0), is 
unambiguous since it is just obtained from following the 
dispersion of the EDC peaks which clearly fall on the 
umklapp band dispersion in Fig. 12(c). 

The Fermi crossing U4 requires more care since the 
umklapp and main band dispersions are almost degen- 
erate in the vicinity of (tt, 0). To disentangle these two 
contributions, we exploit the polarization selection rules 
discussed in Section IV. In Fig. 12(c) we use filled cir- 
cle symbols to denote data obtained in an odd geometry, 
i.e., initial state odd under reflection in the corresponding 
mirror plane, and open circles to denote even geometry. 
We see from the dispersion plotted in Fig. 12(c) that the 
main band signal is seen in the even geometry, since it is 
a dipole-allowed transition. (Actually in this polarization 
both the main and the umklapp bands should contribute 
but in the 17-22 eV photon energy range the main 
band intensity is much larger than that of the umklapps.) 
However in the odd polarization (filled circles along (0, 0) 
to (tt, 0)) the main band is dipole-forbidden and thus the 
weaker umklapp band, which does not have any symme- 
try restrictions here, dominates in this geometry. From 
the dispersion, and in particular the polarization geome- 
try in which it is observed, we clearly see that the Fermi 
crossing U4 must be associated with the umklapp Fermi 
surface, and not with the main band. 

Fig. 13 shows various cuts at a photon energy of 22 eV 
for Bi2212-OD87K. The Fermi crossings are determined 
using the symmetrization method and the Fermi surface 
is found to be a hole barrel centered at (7r,7r). Notice 
that in each cut umklapp bands (broken lines) can be 
identified. The labels Ul to U3 correspond to particular 
Fermi crossings of the umklapp band as shown in Fig. 12. 
In particular along the (0, 0) to (tt, tt) direction (right-side 
panels) the Ul and U2 crossings of the umklapp bands 
can be clearly observed. 

These umklapp bands are also responsible for the 
Lu — peaks observed in the symmetrized data beyond 



Icf in the other panels. The new data, obtained at much 
higher density in the zone, allows us to directly visu- 
alize the main Fermi surface, together with the ghost 
Fermi surfaces due to umklapp bands, using the same 
procedure as in Fig. 11. The magnitude of logarithmic 
gradient |Vk/(k)|/J(k) for Bi2212-OPT90K over part of 
the Brillouin zone in the Y quadrant (defined in Fig. 12) 
is shown in Fig. 14. From the intensity pattern in this 
plot, one can clearly see the main Fermi surface in the 
middle, which is a large hole-like barrel, and also one of 
the umklapps (U3 in Fig. 12 notation). The other umk- 
lapps are weaker in intensity, but would be visible in the 
figure if a log intensity scale had been used. It is quite 
satisfying to see indications of all the features deduced 
earlier (using other methods) in the plots obtained using 
the straightforward logarithmic gradient method. 

The analysis of the photon energy dependence of the 
ARPES data in Bi2212 shows a similar trend to the one 
observed in Bi2201. As an example of this, we plot 
in Fig. 15 the photon energy dependence of the ratio 
/(O.Ttt, 0)//(7r, 0) of the energy integrated intensity mea- 
sured at (O.Ttt, 0) and (tt, 0) for OD88K. Both these k 
points are inside the occupied part of the zone, and we 
would not expect the momentum distribution n(k) to 
vary significantly from one point to the other. Thus, 
following Eq. (jg), any significant deviation of this ratio 
from unity as a function of incident photon energy must 
be attributed to the matrix elements. While around 20 
eV, /(O.Ttt, 0)//(7r, 0) is close to unity indicating a small 
k-dependence of the matrix elements in this part of the 
zone, the ratio peaks to about 2.5 at 38 eV, signaling the 
suppression of intensity around (tt, 0), similar to what is 
seen in Bi2201. Even more interesting is the observation 
that around 54 eV the /(O.Ttt, 0)//(7r, 0) becomes much 
smaller than one. Fig. 15 illustrates once more how dan- 
gerous it would be to infer a Fermi crossing from intensity 
variations alone. ■— ■ 

We have further observedll3 that at photon energies 
close to 30 eV, where the main band is strongly sup- 
pressed around (tt, 0), the superlattice contributions are 
in fact strongly enhanced in this region. The reason is 
that these superlattice intensities originate from regions 
of reciprocal space where the matrix elements are less 
suppressed. 

To summarize, it is very important that the superlat- 
tice contributions be differentiated from the main band 
using polarization selection rules, and all Fermi crossings 
carefully checked by a combination of symmetrization 
analysis together with careful studies of the integrated 
intensities as shown above. 



X. CONCLUSIONS 



high Tc copper oxide based materials, Bi2201 and Bi2212. 
However we believe that these methods should prove to 
be useful for a large class of quasi 2D materials. 

The high Tc materials are hard to analyze because of 
the absence of sharp quasiparticle peaks in their nor- 
mal state, and anomalously weak dispersion in parts of 
the zone. Nevertheless, the symmetrization method dis- 
cussed in this paper is able to deal with both these issues. 
It effectively removes the Fermi function from the EDC 
and determines kj^ as that point in k-space at which the 
spectral function peaks at the chemical potential. 

It is very useful to supplement this analysis with stud- 
ies of the momentum distribution, however one has to be 
very careful about matrix element effects. We show that 
by analyzing data at different incident photon energies, 
one can unambiguously distinguish between loss of inte- 
grated intensity arising from matrix element variations 
from that due to genuine structure in n(k). In this con- 
nection we have also shown the usefulness of studying the 
gradient of the logarithm of the integrated intensity. 

We emphasize that not recognizing the role of matrix 
elements can lead to paradoxical conclusions like changes 
in Fermi surface topology with photon energy, which of 
course makes no sense. Some authors have ascribed the 
differences between the 22 eV and 34pey_|photon energy 
data to additional states around (tt, 0)E3li3. We find that 
there is no necessity to invoke such states. 

Bi2212 has an additional complication arising from the 
effect of the BiO superlattice modulation on the ARPES 
data. However, as we argue above, the Fermi surface 
crossings arising from "ghost" images (superlattice umk- 
lapp bands) can be clearly differentiated from the intrin- 
sic planar Cu02 Fermi surface by exploiting polarization 
selection rules. 

In conclusion, all our data on Bi2201 and Bi2212 when 
analyzed using the methods described above indicate a 
single large Fermi surface centered around (tt, tt) inde- 
pendent of the photon energy. 
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In this paper we have carefully enunciated the criteria 
to be used in determining the Fermi surface from ARPES 
data. We have illustrated these ideas using data on two 



APPENDIX A: SYMMETRIZATION AND 
RESOLUTION 

In this Appendix we discuss in detail the effect of 
experimental resolution on the symmetrization method. 
For simplicity, we discuss the elimination of the Fermi 
function from the EDC at hp, which requires an assump- 
tion of particle-hole symmetry. As discussed in Section 
V, no such assumption was required for the determina- 
tion of \ip via symmetrization. 

Symmetrization was first introduced by us in Ref. Eq 
and used extensively for studying the self energy in 
Ref. E% The main result 

/sy,n(kF,t^) = /(kF,Cj) +/(kF,-w) = IoA(kF,U!). 

(Al) 
follows immediately from Eq. (0) by using the identity 



fi-Lu) = 1 - f{u) 



(A2) 



obeyed by the Fermi function, together with the assump- 
tion of particle-hole symmetry at low energies (w less 
than few times the temperature): 



A(\<iF,^) = A(k_F,-w). 



(A3) 



Let us now see how symmetrization works in the pres- 
ence of finite energy and momentum resolutions. For 
clarity of presentation, we discuss these one at a time, 
although both can be trivially treated together. With a 
finite energy resolution, Eq. (0) is generalized to 



/(k,c^) = /o 



+ 00 



dw'R{uj - uj')f{uj')A{\i, w'), (A4) 



where R is typically taken t o be a Gau ssia n. Using R{lu — 
Lo') = R{uj' — u)) and Eqs. (A2) and ( |A3|) , we can easily 
see that 



/sym(k_F, Uj) = /(kf, Uj) + /(k_F, -w) 



/o 



+ OC 



dLu'R{uj-uj')A{kF,Lj'). 



(A5) 



Thus symmetrization succeeds in removing the effect of 
the Fermi function from inside the convolution integral. 
Next consider the effect of a small, but finite, k- 
window. In its presence, Eq. (0) is replaced by 



/(k,c.) = /o/HV'A(k',^). 



(A6) 



where ^ is shorthand for summation over k' within the 
window centered about k. We ignore the k-variation of 
the prefactor Iq within this small window which allows 
us to pull it out of the sum. Next we need to extend our 
particle-hole symmetry assumption to k's slightly away 
from hp. We require 



A(ek,a;) = yl(-ek, -w) 



(A7) 



valid for |a;| and |ek| both less than few tens of meV. 
Note that we have rewritten the first argument of the 
spectral function as ek, which can be linearized in the 
vicinity of hp as ek — vp-Sk. where (5k = (k — k-p). 
The symmetrized intensity is thus given by IgymO^F , '-^) — 
lofiio) E' M^F-Sk, ij) + Iof{~io) E' A{wp.Sk, -iv). We 
assume a symmetric window, so that if kp + dk is within 
the window, then so is kp — 6k. Thus we can rewrite the 
second term as E Ajvp .Sk, —uj) = E A{—vp.6k,—u!), 
which on using E q. (A7 ) is given by E A{-v p .dk, lo) . Fi- 
nally, using Eq. ([A^), we get 



,{kp,uj)^ioy 'A{k',io) 






(A8) 



Combining the arguments that led to Eqs. (A5) and 
( |A§| ), we see that the symmetrization procedure works in 
the presence of bothr-aiiergy and momentum resolution. 
In our previous workEa-EJ we had used this procedure to 
analyze data at kp mostly in the pseudogap and super- 
conducting states. 
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FIG. 1. Rocking curve of a (0,0,10) reflection on a Bi2212 
sample showing the large structural coherence length. 



FIG. 3. (a) Electronic structure of the valence band of 
Bi2212 obtained by taking the second derivative of EDCs such 
as those shown in (b) at hu =22 eV. 



FIG. 4. (a) Arrangement of the photon beam and detector 
in order to make use of the photoemission selection rules, 
(b) Parity of the Cud^2_y2 orbitals hybridized with the 02p 
orbitals. (c) EDCs showing the parity of the orbitals in (b) 
obtained at hiy =22 eV. 



FIG. 5. Determination of the Fermi crossing along 

(0,0) -^ (7r,7r) in Bi2212-OD88K obtained at hu =22 eV, 
T=180 K: (a) EDCs with the state at the Fermi momentum 
shown in bold; (b) Symmetrization of the EDCs in (a); (c) 
Dispersion obtained from the data in (a). 



FIG. 6. EDCs from Bi2201-OD23K obtained at hu =22 
eV, T=25 K along cuts perpendicular to the (0, 0) -^ (tt, 0) 
direction, with kx indicated in each panel. 



FIG. 7. (a,b,c) Symmetrization of selected data in Fig. 6. 
(d) EDCs along the (0, 0) — > (27r, 0) direction; (e) symmetriza- 
tion of the EDCs in (d), showing no Fermi surface crossing. 



FIG. 8. Determination of the Fermi crossing in 

Bi2201-OD23K at hu =34 eV, T=25 K. (a) EDCs along the 
(0, 0) — > (tt, 0) -^ (tt, tt) directions. Symmetrized data (b) 
along (0, 0) ^ (tt, 0) showing no Fermi crossing and (c) along 
(7r,0) -^ (tt, tt) showing a clear Fermi crossing. 



FIG. 9. Bi2201-OD23K, (a) Integrated intensity at 

hu =22 and 34 eV along the (0,0) -^ (7r,0) -^ (7r,7r) di- 
rections, (b) Comparison of the ARPES lineshape measured 
at 22 (dashed lines) and 34 eV (solid lines) at three different 
k points. 



FIG. 10. (a) Model calculation of the dispersive band 
along the (0, 0) — > (tt, 0) direction, assuming no (dashed lines) 
and strong (solid lines) k-dependence of the matrix elements. 
The variation of the matrix element observed in Fig. 9 at 
34 eV is simulated by a sin'*(0.6fca;) function, (b) Integrated 
intensity over a narrow (-50,-f 50 meV) and wide (-350, +50 
meV) energy range using the model calulation shown in (a), 
(c) Integrated intensity over the same integration ranges as 
in (b) obtained experimentally on Bi2201-OD23K at hv =34 
eV. 



FIG. 2. Energy distribution curves (EDCs) obtained at 
hv =22 eV for a Bi2201-OD4K sample, showing the symmetry 
about the (0,0) ^ (27r, 0) (a) and the (tt, -tt) ^ (7r,7r) (b) 
directions. 
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FIG. 11. Bi2201-OD0K: (a and c) Integrated intensity 
(over -310 to +90 meV) /(k) measured at hu =22 and 
28 eV, T=16 K around the (tt.O) point. Notice that the 
intensity maximum depends strongly upon the photon en- 
ergy hu. (b and d) Corresponding gradient of the loga- 
rithm, |Vk/(k)|/J(k), the maxima (white) which correspond 
to Fermi crossings and clearly show that, independent of the 
photon energy, the Fermi surface consists of a hole barrel cen- 
tered around (7r,7r). 



FIG. 16. Plot of symmetrized intensities obtained from 
simulation for five k points {(p = 2.5° — 6.5°) for a. 6 = 14° 
cut in k-space {TM geometry). The curve closest to kp 
(= 4.6°) is shown in bold. The dispersion was chosen from 
the tight-binding fit of Ref. pA. A (constant) linewidth broad- 
ening of 50 meV was used, with a Gaussian energy resolution 
of (T = 15meV, and a 1° radius k-window. A Fermi function 
with T = 14 K was used, and matrix elements were ignored. 
Using these parameters, the five EDCs were generated and 
then plotted after symmetrization. 



FIG. 12. (a) In the presence of the superstructure in the 
BiO layers, the outgoing electrons can be diffracted, thus giv- 
ing rise to additional umklapp bands as shown in (b and c). 
(b) Main (thick) and umklapp (thin) Fermi surfaces. Selected 
Fermi crossings of the umklapp bands relevant for Figs. 13 
and 14 are labeled from Ul to U5. (c) Dispersions obtained 
from earlier measurements on Bi2212-OD87Ktl. Filled cir- 
cles denote data obtained in an odd polarization, i.e., initial 
state odd under refiection in the corresponding mirror plane. 
Open circles denote even polarization, and open triangles cor- 
respond to a mixed polarization. 



FIG. 13. Bi2212-OD87K. The top panels show the EDC 
cuts taken at hi/ =22 eV, T=100 K together with the po- 
larization geometry used. Various EDCs are shown in the 
middle panels, together with the corresponding symmetrized 
data in the lower panels. The curves corresponding to the 
Fermi crossing of the main and umklapp bands are shown 
with thick and broken lines, respectively. The labels Ul to 
U3 correspond to particular Fermi crossings of the umklapp 
band as shown in Fig. 12. 



FIG. 14. Gradient of the logarithm of the integrated inten- 
sity (over -320 to +80 meV), |Vk/(k)|//(k), around (7r,0) for 
Bi2212-OPT90K taken at /ii/=22 eV. Note the large hole-like 
Fermi surface corresponding to the main band. In addition, 
the umklapp band (U3 in Fig. 12 notation) can be seen as 
well. The other umklapps are weaker in intensity, but would 
be visible if a log intensity scale had been used instead. 



FIG. 15. Photon energy dependence of the integrated in- 
tensity ratio (over -600 to +200 meV) /(0.77r,0)//(7r, 0) of 
Bi2212-OD88K. The ratio is near unity at 20 eV, increases at 
30 eV to two and larger, and decreases at 54 eV to less than 
one-half. This figure illustrates how dangerous it would be to 
infer a Fermi crossing from intensity variations only. 
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